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Abstract 

The chiral-odd transversity generalized parton distributions (GPDs) of the nucleon can be 
accessed experimentally through the exclusive photoproduction process "/ + N — !> n + p + N' , in 
the kinematics where the meson pair has a large invariant mass and the final nucleon has a small 
transverse momentum, provided the vector meson is produced in a transversally polarized state. 
We calculate perturbatively the scattering amplitude at leading order in a^. We build a simple 
model for the dominant transversity GPD Ht{x, ^, t) based on the concept of double distribution. 
We estimate the unpolarized differential cross section for this process in the kinematics of the 
Jlab and COMPASS experiments. Counting rates show that the experiment looks feasible with 
the real photon beam characteristics expected at JLab@12 GeV, and with the quasi real photon 
beam in the COMPASS experiment. 



1 Introduction 



Transversity quark distributions in the nucleon remain among the most unknown leading twist hadronic 
observables. This is mostly due to their chiral odd character which enforces their decoupling in most 
hard amplitudes. After the pioneering studies [T], much work [2] has been devoted to the exploration 
of many channels but experimental difficulties have challenged the most promising ones. 

On the other hand, tremendous progress has been recently witnessed on the QCD description 
of hard exclusive processes, in terms of generalized parton distributions (GPDs) describing the 3- 
dimensional content of hadrons. Numerous experimental and theoretical reviews [3] exist now on this 
quickly developing subject. It is not an overstatement to stress that this activity is very likely to 
shed light on the confinement dynamics of QCD through the detailed understanding of the quark and 
gluon structure of hadrons. 

Access to the chiral-odd transversity generalized parton distributions 4 , noted Ht, Et, Ht, Et, 
has however turned out to be even more challenging |5j than the usual transversity distributions : one 
photon or one meson electroproduction leading twist amplitudes are insensitive to transversity GPDs. 
A possible way out is to consider higher twist contributions to these amplitudes [B] , which however are 
beyond the factorization proofs and often plagued with end-point singularities. The strategy which we 
follow here, as initiated in Ref. [71 [S], is to study the leading twist contribution to processes where more 
mesons are present in the final state; the hard scale which allows to probe the short distance structure 
of the nucleon is now the invariant mass of the meson pair, related to the large transverse momentum 
transmitted to each final meson. In the example developed previously [T] [5] , the process under study 
was the high energy photo (or electro) diffractive production of two vector mesons, the hard probe 
being the virtual "Pomeron" exchange (and the hard scale being the virtuality of this pomeron), in 
analogy with the virtual photon exchange occuring in the deep inelastic electroproduction of a meson. 
A similar strategy has also been advocated recently in Ref. [S] to enlarge the number of processes 
which could be used to extract information on chiral-even GPDs. 

The process we study here ([TP]) 

-f + N^Tr++p°T + N', (1) 

is a priori sensitive to chiral-odd GPDs because of the chiral-odd character of the leading twist 
distribution amplitude of the transversally polarized p meson. Its detailed study should not present 
major difficulties to modern detectors such as those developed for the 12 GeV upgrade of Jlab or for the 
Compass experiment at CERN. The estimated rate depends of course much on the magnitude of the 
chiral-odd generalized parton distributions. Not much is known about them, but model calculations 
have been developed in [5] for the ERBL part and in Ref. [TTl [T2J [T3]; moreover, a few moments 
have been computed on the lattice [13]. To supplement this and use the recent phenomenological 
knowledge acquired on the transversity quark distributions through single inclusive deep inelastic 
data, we propose in this paper a parametrization of the (dominant) transversity GPD based on 
the concept of double distributions. 

Let us now explain how we factorize the amplitude of this process and what is the rational of this 
extension of the existing factorization proofs in the framework of QCD. The basis of our argument is 
two-folded. 

• We use the now classical proof of the factorization of exclusive scattering at fixed angle and large 
energy |15j . The amplitude for the process 7+7r— >-7r-|-/ois written as the convolution of mesonic 
distribution amplitudes and a hard scattering subprocess amplitude j+{q + q) — > iq + q) + iq + q) 
with the meson states replaced by coUinear quark- ant iquark pairs. This is described in Fig. [TJi. 
The demonstration of absence of any pinch singularities (which is the weak point of the proof 
for the generic case A + B ^ C + D ) has been proven in the case of interest here US] . 

• We extract from the factorization procedure of the deeply virtual Conipton scattering amplitude 
near the forward region the right to replace in Fig. [T^ the lower left meson distribution amplitude 
by a iV — 7^ TV' GPD, and thus get Fig.[T}). Indeed the same coUinear factorization property bases 
the validity of the leading twist approximation which either replaces the meson wave function 
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by its distribution amplitude or the N N' transition to its GPDs. A slight difference is 
that light cone fractions (z, I — z) leaving the DA are positive, but the corresponding fractions 
{x + — x) may be positive or negative in the case of the GPD. The calculation will show that 
this difference does not ruin the factorization property, at least at the order that we are working 
here. 

One may adopt another point of view based on an analogy with the timelike Compton scattering 

-fN ^-/*N' ^ ^i+^i^N' , (2) 

where the lepton pair has a large squared invariant mass Q^, is instructive. This process has been 
thoroughly discussed [T^ in the framework of the factorization of GPDs, and it has been proven that 
its amplitude was quite similar to the deeply virtual Compton scattering one, being dominated at 
lowest order by the handbag diagram amplitude convoluted with generalized quark distributions in the 
nucleon. There is no ambiguity in this case for the definition of the hard scale, the photon virtuality 
Q being the only scale present. Although the meson pair in the process ([T]) has a more complex 
momentum flow, we feel justified to draw on this analogy to ascribe the role of the hard scale to the 
meson pair invariant squared mass. However, to describe the final state mesons by their distribution 
amplitudes (DAs), one needs in addition a large transverse momentum (and thus a large Mandelstam 
t', see Fig.lT}). Practically, we consider kinematics in which \u'\ ^ \t'\ ^ ~ M^^ = (p^ +pp)^. We 
cannot prove, at the level of our study, that Af^^ is the most adequate hard scale. Indeed, applying 
a definite strategy to define a factorization scale requires at least a next to leading (in the strong 
coupling) analysis [18j and this is clearly a major work to be undertaken. 

For both point of view, in order for the factorization of a partonic amplitude to be valid, and the 
leading twist calculation to be sufficient, one should avoid the dangerous kinematical regions where 
a small momentum transfer is exchanged in the upper blob, namely small t' = (p^r — P-/)'^ or small 
u' — {pp — P'y)'^, and the regions where strong interactions between two hadrons in the final state 
are non-perturbative, namely where the invariant masses, Af^^, = (p^r +Pn')'^, ^pN' ~ (Pp +P-/V')^ 
and M^p, are not large enough to suppress final state interactions. We will discuss the necessary 
minimal cuts to be applied to data before any attempt to extract the chiral odd GPDs. However, 
although the ultimate proof of the validity of the factorization scheme proposed in this paper is based 
on comparison of the predictions with experimental data, on the theoretical side it requires to go 
beyond Born approximation considered here which is beyond the scope of the present work. 

Our paper is organized as follows. In section [51 we clarify the kinematics we are interested in and 
set our conventions. Then, in section |3l we describe the scattering amplitude of the process under 
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study in the framework of QCD factorization. Section S] is devoted to the presentation of our model 
chiral-odd GPDs. Section [5] presents our results for the unpolarized differential cross section in the 
kinematics of two specific experiments : quasi-real photon beams at JLab where S'^at ~ 14-20 GeV^ 
and Compass at CERN where S^n ^ 200 GeV^ 

As a final remark in this introduction, let us stress that our discussion applies as well to the 
case of electroproduction where a moderate virtuality of the initial photon may help to access the 
perturbative domain with a lower value of the hard scale Mj^p. 

2 Kinematics 

We study the exclusive photoproduction of a transversely polarized vector meson and a pion on a 
polarized or unpolarized proton target 

7(g) + N{pi,\) -> 7r{p^) + pt{Pp) + N'[p2, A') , (3) 

in the kinematical regime of large invariant mass M^p of the final meson pair and small momentum 
transfer t = (pi — ^2)^ between the initial and the final nucleons. Roughly speaking, these kinematics 
mean a moderate to large, and approximately opposite, transverse momentum of each meson. Our 
conventions are the following. We decompose momenta on a Sudakov basis as 

A;^==an^ + 6p^ + /s^, (4) 

with p and n the light-cone vectors 

p'^ = ^(1,0,0,1) n^=^^(l, 0,0,-1) P-n^'-, (5) 

and 

A:'^ (0,r,F,0), ki = -k1. (6) 

The particle momenta read 

p; = (l+0p" + j7^^n" , = (1 - P" + t f )' + ^1 . «" = "". (7) 

as 2 

^ apS 2 

with a — I — a and M, m^r, rup the masses of the nucleon, the pion and the p meson. From these 
kinematical relations it follows 

r, , ^ jPt - ^At)^ + ml ^ {pt + ^At)^ + ml 
s a sap 

and 



l-a-ap = -j^-:^^+-j^^^. (10) 



2^NP A'f 

The total center-of-mass energy squared of the 7-N system is 

S^N^{q+Pif ^il + Os + M\ (11) 

^ is the skewedness parameter which can be written in terms of the r variable used in lepton pair 
production, as 
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On the nucleon side, the transferred squared momentum is 

t=(P2-Pi)' = -Y3|A2-^S__. (13) 

The other various Mandelstam invariants read 



s' ^ (P.+Ppf - Mj, = 2e.(l-^^^)-A?i±|, (14) 

a 

u = -{pp-q) = (16) 



and 





s(i-0/ V 2 



(18) 

The hard scale M^^ is the invariant squared mass of the (tt^, system. The leading twist 
calculation of the hard part only involves the approximated kinematics in the generalized Bjorken 
limit: neglecting in front of p± as well as hadronic masses, it amounts to 

Ml « 4, (19) 
ap ~ 1 — a = a , (20) 

r - a ^ (21) 



M2 



-t' « a Mlp and - u' « a M^^ . (22) 

The typical cuts that one should apply are —t',—u' > and M^j^, — {p.^ + pn')'^ > M^, 
MpN> — [Pp + Pn'Y > where A ^ ^qcd and Mp^ is a typical baryonic resonance mass. This 
amounts to cuts in a and a at fixed M^^, which can be translated in terms of u' at fixed and t. 
These conditions boil down to a safe kinematical domain {—u')min < —u' < {—u')max which we will 
discuss in more details in Section [5] 

In the following, we will choose as kinematical independent variables . 



3 The Scattering Amplitude 

We now concentrate on the specific process 

7(q) +p(pi. A) ^ 7r+(p,) + p"r{Pp) + n{p2, A') . (23) 

Let us start by recalling the non-perturbative quantities which enter the scattering amplitude of our 
process (123^ . The transversity generalized parton distribution of a parton q (here q — d) in the 
nucleon target at zero momentum transfer is defined by jl^ 

< n{p2,X')\ a+^Y'u (I) |p(pi. A) > 

= w(p2,A')CT+^7^u(pi,A)y dxe-i'^'-P^+Pi^y Hf{x,^,t), (24) 
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where A and A' are the Hght-cone hehcities of the nucleons p and n. Here Hlj^'^ is the flavor non-diagonal 
GPD [T^ which can be expressed in terms of diagonal ones as 

= H^-H^. (25) 

The chiral-odd light-cone DA for the transversely polarized meson vector pj^, is defined, in leading 
twist 2, by the matrix element [20] 







where t^(j)p) is the p-meson transverse polarization and with = 160 MeV. 

The light -cone DA for the pion tt^ is defined, in leading twist 2, by the matrix element (see for 
example [31]) 

(0| J(z)7Vw(-^)k+(p)> = du e-*(2«-i)p-. ^^(^)^ (27) 



"'0 

with — 131 MeV. In our calculations, we use the asymptotic form of these DAs : (prriu) ~ 4>±{u) — 
6 uu. 

We now pass to the computation of the scattering amplitude of the process (1^ . As the order 
of magnitude of the hard scale is greater than 1 GeV^, it is possible to study it in the framework 
of QCD factorization, where the invariant squared mass of the (7r+, p^) system M^^ is taken as the 
factorization scale. 

The amplitude gets contributions from each of the four twist 2 chiral-odd GPDs Et , Ht , Et , Ht- 
However, all of them but Ht are accompanied by kinematical factors which vanish at Af = . The 
contribution proportional to Ht is thus dominant in the small t domain which we are interested in. 
We will thus restrict our study to this contribution, so that the whole dependence will come from 
the t-dependence of Ht, as we model in Sec. S] Note that within the collinear framework, the hard 
part is computed with Aj = 0. 

Thus we write the scattering amplitude of the process (j23p in the factorized form : 

A{t,Mlp,pT) = ^j ^dx dv dz{T^{x,v,z)-T''{x,v,z))H^\x,£,,t)^^{z)^^{v), (28) 

where T" and T'^ are the hard parts of the amplitude where the photon couples respectively to a u- 
quark (Fig. [2^) and to a d-quark (Fig.[2)D). This decomposition, with the ^ prefactor, takes already 

into account that the p"-meson is described as ^^^'^ ■ 

For this process, one has two kinds of Feynman diagrams : some without (Fig. [5]) and some with 
a 3-gluon vertex (Fig. |3|). In both cases, an interesting symmetry allows to deduce the contribution 
of some diagrams from other ones. This is examplified in Fig. [5J The transformation rules 

X — > —x u ^ u V ^ V e„ — > (29) 

relate the hard amplitude of Fig. [^b to the one of Fig. This reduces our task to the calculation of 
half the 62 diagrams involved in the process. 

Let us sketch the main steps of the calculation on the specific example of the diagram of Fig. [2^, in 
the Feynman gauge. Using the notation ^ — k^"f''\ the amplitude reads : 

T2a{x,v, z) = Tr[{if^i)^-i^){-ig-f'')f{p'2 + vpp + zp^){ieui{q))f{p'i - ^Pp - zp^) 

X TrcrM]j^^G{p',+vpp)G{vpp + zp^), (30) 
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Figure 3: Representative diagram with a 3 gluon vertex. 



where the fermion propagator is (we put all quark masses to zero) : 
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and 



fc2 



is the gluonic propagator. Trc is the trace over color indices and the factors 
from Fierz decompositions. The corresponding expression for the diagram [^b 



and 



4Afc 



(32) 



come 



T2b{x,v,z) 



Tr[(*/.M')(-*.97^)(2*^">^/p^)(-*.97'')(^"''7') 
{-i91i^)f{p'2 + vpp + zp^){ied^{q))f(p'-^ ~ vpp - zp^)(-z.g7^)] 

■G{~p[ + vpp)G{-vpp - zp^) 



(8iv,)2 



{Nt ■ sum ■ ?^t) - {Nt ■ Pt){e^t ■ a^pt) + ^^{N^ ■ e,t){pt ■ a;,) 



zvv[{az + av){x + ^ — «e) — 2£^zv\ 



(33) 



justifies the symmetry we quoted a few hnes above. Thus the hard part of the diagram [2^ is propor- 
tional to 



T. 



1 



2a 



[(^2 + vi)p + zi)^)"^ + ie] [(/5'i - vi)p - zi)^Y + ie] [{p'^ + vppY + ie] [(vpp + zp^) 



(34) 



and the ie prescription in the 4 propagators leads to the fact that the scattering amplitude gets both a 
real and an imaginary parts. Integrations over v and z have been done analytically whereas numerical 
methods are used for the integration over x. The first integration is rather straightforward. The 
second integration is more involved because of the presence of ie terms inside the integrand, and in 
particular as an argument of logarithmic funtion, leading in the final result to appearance of imaginary 
parts. For example, the integration over z of T^j^ requires to evaluate integrals of the type 



dz — — 

2iz-aX 



loa 



aXz 



aX + z {aX - 20 



(35) 



where X = x — ^ + ie contains all the dependence of the integrand on ie. 

Nevertheless, since we have rewritten the x-dependence of propagators with the new variable X, it is 
possible to calculate this integral analytically without any problem. Thus the expression (1551) gives 



1 Ll2 

\2i 2i 



1 - 



aX 



1 



aX 



Li? 

2e 



1 



aX 



Lio 

2^ 



1 



aX 



(36) 



Lorentz invariance and the linearity of the amplitude with respect to the polarization vectors and 
with respect to the nucleons' spinors allow us to write the amplitude as : 



A = (elipp) ■ N^,,J{e^i_ ■ Pt)A' + (eUPp) ■ e^i-)iN^,x. ' Pt)B' 

+ (elipp) -prKNt,^, ■ e,^)C' + (el(Pp) • PT)(A^i,A. ■ PT){e^^ ■ Pt)D' 

+ {elipp) ■p){Ni^^^_ ■ e,i:)E' + [e^ipp) ■ p)iN^^^^ ■ pT){e^^ ■ pt)F\ (37) 

where A' , B' , C , D' , E' , F' are scalar functions of s, ^ a and M^p, and the transverse polarization 
of p-meson 

''^M={^.^^.^,J0^f.) (38) 

is expressed in terms of e± — — -^(±l,i,0). e^^^ is the transverse polarization of the on-shell photon 
and 

- ^9Tu{p2. A2)^7.7'«(Pi, Ai) (39) 
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is the spinor dependent part which expresses the nucleon hehcity flip with g'^ = diag(0, —1, —1, 0). 
To be more precise, the expressions of this 2— dimensional transverse vector read 



= -4*^1-^2(0,1,0,0) A^f|'+£ = 4^1-^2(0, 0,1,0) 



N 



-4v/r^(0, 0,1,0) N_ 



4iv/l -^2(0, 1,0,0) 



(40) 
(41) 



assuming that these nucleons are polarized along the x axis. 

Since the DA of pj. (Eq. (^5]) ) introduces the factor ej.{pp)pp — e'j.{pp)p'^, any term proportional to 
p^ in its polarisation does not contribute to the amplitude. On may then replace 



2a 



2a 



Pt ■ e± 
a'^s + p^ 

Pt ■ e± 

a^s + Pf 



ip' 



1 



(0,e±) 



Pt 



a-^s+Pf 



(0,e±). 



(42) 



Consequently, the amplitude of this process can be simplified as 

A - {Nfel){pfe^t)A+{Nfe^t){Pfel)B 

+ [Nt ■ PtW^t ■ el)C + (Nt ■ Pt){Pt ■ e^tWt ■ Sl)D (43) 
where A, B, C, D are also scalar functions of s, ^, a and M^p. 

The final result for each particular diagram is rather lenghty, and because of that we do not present 
explicit final results for scalar functions A, B, C, D of 



4 Transversity GPD and Double Distribution 



In order to get an estimate of the differential cross section of this process, we need to propose a model 
for the transversity GPD H!^{x,£,,t) {q — u, d). Contrary to what Enberg et al. have done [S], here 
we must get a parametrization in both ERBL (] — f;^[) and DGLAP ([—1; — ^] IJ [f; 1]) a;— domains. 
We use the standard description of GPDs in terms of double distributions [15] 



ir|,(x,C,t = 0) / d(3da5{(i + ^a-x)f^{(i,a,t^Q). 
Jn 



(44) 



where /|. is the quark transversity double distribution and $7 = + |a| ^ 1} is its support domain. 
Moreover we may add a D-term contribution, which is necessary to be completely general while 
fulfilling the polynomiality constraints. Since adding a D-term is quite arbitrary and unconstrained, 
we do not include it in our parametrization. We thus propose a simple model for these GPDs, by 
writting /|. in the form 



/«(/?, a,t = 0)= n(/3, a) S q(/3)e(/3) - n(-/3, a) 5q{-P) , 



(45) 



_ 3 (1- 



3 — is a profile function and 5q^ Sq are the quark and antiquark transversity 



where n(/3, a) - 4— (T3^ 
parton distribution functions (PDF). The transversity GPD H!^ thus reads 



H^{x,^,t^O) = Q{x>0 



1-5 



dy . 
1+^ 4 



3{l-x + ^yf-y^ 



(l^x + ^y)^ 



6q{x ~ ^y) 



+ e(e>x>-e) 



T 3 il-x + ^yy-y\ 
4 (1 — a; + ^yy 



^^^3(1+0;- ^yf 



4 (l+.-Cy)3-'^^"^-" + ^^) 



- ei-^>x) 



1+e 



l + j: 

1-e 



dy 



3il + x-^yf 



4 {1 + x-^y)^ 



-Sq{~x + ^y). 



(46) 
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For the transversity PDFs Sq and Sq, we use the parametrization proposed by Ansehnino et al. [23] 



5u{x) = 7.b-0.b-{l-xf ■{xu{x)^x^u{x)), (47) 

5u{x) = 7.b-0.b-{l-xf ■(xu{x)^xl\u{x)), (48) 

5d{x) = l.b-{-Q.&)-{l-xf -{xdix^+xMix)), (49) 

8d{x) = l.b-{-Q.Q)-{l-xf ■{xd{x)+xM{x)), (50) 

where the hehcity-dependent PDFs l^q{x), l^q{x) are parametrized with the help of the unpolarized 
PDFs q{x) and q{x) by [H] 

AM(a;) = y/x ■ u{x) , (51) 

Am(x) = -0.3 • 2;°-'' • u(a;) , (52) 

M{x) = -^0.7 ■ y/^- d{x) , (53) 

Ad{x) = -0.3 ■ x^ "^ ■ d{x) , (54) 



and the PDFs q{x), q{x) come from GRV parametrizations [25]. Ah these PDFs are calculated at 
the energy scale /i^ = 10 GeV^. Fig. |3] represents H^{x,(,,t = 0) and H!f(x,£,,t = 0), respectively, 
for different values of which are determined through ([T^ for S^n = 20 GeV^ of JLab and for M^^ 
equal 2, 4, 6 GeV^. Similarly, Fig. [5] represents H^{x,£_,t — 0) and H^{x,£,,t = 0), respectively, for 
different values of ^, which are determined through ([T^ for S'-^at = 200 GeV^ of Compass and for Af^^ 
equal 2, 4, 6 GeV^. 

These two GPDs show some common features like a peak when x is near their order of magnitude 
and the fact that they both tend to zero when x tends to ±1. The main difference is their opposite 
sign. The restricted analysis of Ref. [8] based on a meson exchange is insufficient for this study since 
it only gives us the transversity GPDs in the ERBL region. The MIT bag model inspired method 
of Ref. [iT] underestimates the value of Ht[x,£) in the ERBL domain because this model does not 
take into account antiquark degrees of freedom. One can notice that these GPDs have the same order 
of magnitude but some differences with other models like light- front constituent quark models |12| . 
principally due to the fact that in parametrizations have been done at 0.1 GeV^ whereas 

our model is calculated at /i^ ~ 10 GeV^. Other model-studies have been developed in the chiral 
quark soliton model and a QED-based overlap representation |13) . 

The t-dependence of these chiral-odd GPDs - and its Fourier transform in terms of the transverse 
localization of quarks in the proton [2 6) - is very interesting but completely unknown. We will describe 
it in a simplistic way as: 

H;j.{x,^,t)^ H^ix,^,t = 0) X Fnit), (55) 

where 

FhH) = ^ (56) 

is a standard dipole form factor with C = .71 GeV^. Let us stress that we have no phenomenological 
control of this assumption, since the tensor form factor of the nucleon (in fact even the tensor charge) 
has never been measured. 
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(a) (&) 

Figure 4: Transversity GPD H^{x,^,t = 0) (a) and H^{x,i,t = 0) (b) of the nucleon for C = -HI 
(solid line), ^ = .176 (dotted line), ^ — .25 (dashed line), corresponding respectively to M'^p/S^M 
equal to 4/20, 6/20 and 8/20. 




(«) (&) 

Figure 5: Transversity GPD H^{x,£„t = 0) (a) and H^{x,^,t = 0) (b) of the nucleon for f = .01 
(sohd line), ^ = .015 (dotted line), ^ = .02 (dashed line), corresponding respectively to Al^p/S^jy 
equal to 4/200, 6/200 and 8/200, the lower plots being blow-ups of the central part of the upper ones. 
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5 Unpolarized Differential Cross Section and Rate Estimates 




(a) 



-u'(GcV2) 



-u'(GeV2 



Figure 6: Variation of the differential cross section (1551) (nb.GeV ^) with respect to \u'\ at M^^ = 
3 GeV^ (a) and M^^ = 6 GeV^ (b) with S^n = 20 GeV^. The hnes on the left correspond to the 
constraints — u' > 1 GeV^ and A/^^, > 2 GeV^ and the lines on the right correspond to the constraints 
-t' > 1 GeV^ and M^j^, > 2 GcV^ (dashed hne for t = t^in and solid line for t = -0.5 GeV^). 



Starting with the expression of the scattering amplitude (j43p we now calculate the amplitude 
squared for the unpolarized process 



\M\' 



(57) 



Ai A2 



It can seem odd to study the chiral-odd quark content of the nucleon by calculating the cross section 
of an unpolarized scattering but it is enough for now in order to reach this unknown structure. Of 
course it is possible to consider the polarized one by producing the spin density matrix, which will be 
done in a future work. 

We now present the cross-section as a function of t, , —u' which reads 



da 



dtdu' dMlp 



325^^M2^(27r)3- 



(58) 



20 GeV^ for M^^ 
200 GeV2 



We show, in Fig. [SI the differential cross section (|58p as a function of ~u' at S-^n 
= 3 GeV^ i.e. f = 0.085 and for M^^ = 6 GeV^ i.e. ^ = 0.186 and, in Fig. [3 at S-yN 
respectively for M^^ = 3 GeV^ i.e. ^ = 0.0076 and for M^^ = 6 GeV^ i.e. ^ = 0.015. 

To get an estimate of the total rate of events of interest for our analysis, we first get the M^^ 
dependence of the differential cross section integrated over u' and t, 

da__ fO-S 



di-t) 



da 



dtdu'dMlp 



(59) 



The domain of integration over —u' is deduced from the cuts we discussed at the end of section [2] 

M^jy, > 2 GeV^, and is illustrated in Fig. H] Then we get two limits 



^t',-u' > 1 GeV^ and M^^^,, 
-u' : 



for the domain over 



The cuts over —u' and M^j^, give the minimum value for —u' (- 



1 GeV2 or {-u' 



'rmnires 



(res.) } 

which depends on t, S~^n and M^^. For instance, the lines on the left in Figs. |B][7] represent that 
cut : solid lines for t = -0.5 GeV'^ and dashed lines for t — tmin- One can notice that at high 
energy, the nN' system is outside the baryonic resonance region so that —u'^ii^ is always equal 
to 1 GeV^ and for any value of t. 
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0.5 1.0 1.5 2.0 1 2 3 4 



-u'(GeV2) -u'(GeV2) 
(a) (fe) 

Figure 7: Variation of the differential cross section (1551) (nb.GeV"^) with respect to \u'\ at M^^ = 
3 GeV2 (a) and M^^ = 6 GeV^ (b) with ^^^v = 200 GeV^. The sohd line on the left corresponds 
to the constraints — u' > 1 GeV^ and M^^, > 2 GeV^ for any value of t and the lines on the right 
correspond to the constraints —t' > 1 GeV^ and M^^, > 2 GeV^ (dashed line for t = t,nin and solid 
line for t = -0.5 GeV^). 




t 



Figure 8: The phase space domain of integration in the (— t, — w') variables. The upper limit in —u' is 
given by the constraint -u'{t) < {~u')maxit) = -(-t')mm + M^p - t - ml - (with (-i')mm = 1 
GeV^). The lower limit in —u' is given by —u' > 1 GeV^ and —u'{t) > (— u')mm(res.) (i) where 
{—u')jnin{res.){t) IS obtained from the constraint M^j^, > 2 GeV^. The figure illustrates the case 
S-yN = 20 GeV^ and M^^ = 3 GcV^. 



• The cuts over —t' and M^j^, give the maximum value for —u' {{'~u')„iax) which depends on t, 
SjN and M^p. For instance, the lines on the right in Figs. [S][7] represent that cut : solid lines for 
t = -0.5 GeV^ and dashed lines for t = tmin- It is interesting to stress that, for any value of the 
hard scale and of the energy, the pN' system is always outside the baryonic resonance region. 

Moreover, one notices that {—u')„iin decreases and {—u')„iax increases with M^^ at S-yN fixed and 
then the width of the physical region [{^u')min, {~u')max] grows. 

Thus, in Figs. [51 [TUl and [TTl we show the AI^^ dependence of the differential cross section ([5^ . 
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4 5 6 7 8 

M2^(GeV2) 

Figure 9: M^^ dependence of the differential cross section ([5^ (nb.GeV"^) at S-yN = 20 GeV^. 




4 5 6 7 8 



Figure 10: M^^ dependence of the differential cross section (|59| (nb.GeV^^) at S-yM — 100 GeV^. 

Let us first focus on the high energy domain, and discuss the specific case of muoproduction 
with the COMPASS experiment at CERN. Integrating differential cross sections on t, u' and M^p, 
with the cuts specified above and M^^ > 3 GeV^, leads to an estimate of the cross sections for the 
photoproduction of a 7r+p^ pair at high energies such as : 

^7iv^^+p° w {S-,N = 100 GeV^) ~ 3 nb <J^n^^+p'^n' {S-,n = 200 GeV^) ~ 0.1 nb. (60) 

The virtuality of the exchange photon plays no crucial role in our process, and the virtual photo- 
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4 5 6 7 8 



Figure 11: dependence of the differential cross section ([CT)) (nb.GcV^^) at S^n = 200 GeV^. 



production cross section is almost Q^-independent if we choose to select events in a sufficiently narrow 
Q^— window (say Q^^^ < < -5 GeV^), which is legitimate since the effective photon flux is strongly 
peaked at very low values of Q^. The quasi real (transverse) photon flux rlp{Q^, u) reads 



V 2Mp 



27rg2z/2 



1 - 2 



1 - 



El AEfJi. 



(61) 



with the flne structure constant a — 1/137 and Ei the lepton energy (in the laboratory frame). 
Consequently, the rate in a photon energy bin [vi, v-^ corresponding to \S^]si — A5', S-^j^ + A5] with 
S^N = 2 P M = (z^i + V2)M and AS* = 2 Ai/M = {v^ - vx)M is 



a{lN ^Itt+p'^N') 



dQ- 



rV2 

/ r^(g^ i/)cr^.^^^+po^,(Q^i/) 



dQ^ 



dv TlriQ^, y). (62) 



For the muoproduction (£'^ — 160 GeV), one gets the following cross section estimates, flrstly for S-^m 
between 50 and 150 GeV^ 



(J^N^n+p'>N'{S',N^lOOGeV^)x f dQ^ f di^ ^^(g^ z.) 

^ Jo. 02 J25 

IQ-^ nb, (63) 



which yields a rate equal to 3 10 ^ Hz with a lepton beam luminosity of 2.5 10^^ cm ^.s 1, and, for 
SjN between 150 and 250 GeV^ 



aifiN^ pTi+p'^N') ~ cr^jv-^^+pO Ar,(S'7jv = 200 GeV^) x 
~ 5 10"'' nb. 



dg2 



0.02 



125 



75 



dvYt^iQ"^, v) 



(64) 
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which yields a rate equal to 1.3 10~^ Hz with the same lepton beam luminosity. This looks sufficient 
to get an estimate of the transversity GPDs in the region of small ^ of the order 0.01. 

Let us now turn to the lower energy domain, which will be studied in details at JLab. With the 
cuts discussed above and Af^^ > 3 GeV^ , estimates of the cross sections for the photoproduction of a 
■""^Pt P^i^ JLab energies are: 

a^Pf^T,+pO^Pf,{S^N = GeV"^) ~ 15 nh a^]^^^+p0^pf,{SjN = "^0 GeV^) ~ 33 nh. (65) 

In electroproduction {Ee = 11 GeV), applying Eqs. (|6T|) and ((62)) . one gets the total cross section 

a{e-N -> e-TT+p'^N') ~ 0.1 nb. (66) 

Tagging the photons is however required if one aims at a detailed understanding of the reaction and 
at an extraction of the GPD. This is possible at JLab and indeed is well documented for the future 
12 GeV energy upgrade in [37]. More specifically, Hall D will be equipped with a crystal radiator, 
which through the technique of coherent brehmsstrahlung will produce an intense photon beam of 8 
- 9 GeV with an excellent degree of polarization. This leads to the following rate 

~ 5 Hz (67) 

where ^ 10^ photons/s is the photon flux for Hall D and — 1.27 is the number of protons 
per surface in the target (liquid hydrogen of 30 cm), assuming that the efficiency of the detector is at 
100%. 

With a different technique, CLAS12 in Hall B may be equipped with a photon tagger allowing an 
intense (sa 5 10^ photons/s) flux of photons with energy 7 - 10.5 GeV. This will lead to slightly lower 
but still large enough rates. 

Thanks to the high electron beam luminosity expected at JLab, a detailed analysis is possible. 

Moreover, one can make an additionnal comment about the use of the Compass experiment with 
kinematics of JLab, i.e. with photons at low energies. In this context, one gets the following cross 
section estimate for muoproduction for S^n between 20 and 50 GeV^ 

1-1 1-25 

a{^lN ^^^T+ p°,N') ~ a^N^^+p" N'iS-yN = 35 GeV'^) x dQ^ di^ Tt^{Q^ , ly) 

^ Jo. 02 JlO 

~ 0.2 nb , 

which leads to the conclusion that muoproduction with low energy, at Compass, gives greater rates 
(5 10^^ Hz) than with photon with high energy. 



6 Conclusions 

In this paper, we have advocated that the exclusive photoproduction of a meson pair with a large 
invariant mass gives access to generalized parton distributions through the factorization of a hard 
subprocess, provided all the kinematical invariants {s',t',u') which characterize this subprocess are 
large enough. We applied this strategy to access the chiral-odd generalized quark distributions from the 
photoproduction of a tt~^ meson pair with a large invariant mass. We modeled the dominant chiral- 
odd GPD H^'^{x,£^,t) though a double distribution constrained by the phenomenological knowledge 
of the transversity quark distribution hi{x) and hf{x). The calculation of the hard part at the leading 
order in the strong coupling as shows that no divergence nor end-point singularity plagues the validity 
of our approach. From our results, we conclude that the experimental search is promissing, both at 
low real or almost real photon energies within the JLab@12GeV upgraded facility, with the nominal 
effective luminosity generally expected (£ ~ 10'^^ cm^.s"^) and at higher photon energies with the 
Compass experiment at CERN. These two energy ranges should give complementary information on 
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the chiral-odd GPD Ht{x, f , t). Namely, the large ^ region may be scrutinized at JLab and the smaller 
^ region may be studied at COMPASS. 

It is obvious that our model for the chiral-odd GPD may be improved and refined in many 
ways, for instance by adding a D-term which gives a complete parametrization by double distribution. 
Our primary goal in this paper was to prove the feasibility of the study of this physics with this 
physical process. We believe that this task is achieved. 

The described framework opens the way to future studies. Firstly, the contributions proportionnal 
to the other three chiral-odd GPDs should be included in the calculation of the amplitude. Alhough 
they are suppressed by kinematical factors at small t, they constitute an interesting addition to the 
transversity structure of the nucleon. Secondly, it will be interesting to study the polarized process 
and calculate more observables like beam or target spin asymmetries, or like the spin density matrix 
of the vector meson. These two extensions will be discussed in a separate paper. Other two meson 
channels may also be interesting; they deserve a thorough study. 
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